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Abstract 

We apply the method of Bozzolo* Ferrante and Smith for the calculation of alloy energies 
for bcc elements. The heat of formation of several alloys is computed with the help of 
the Connolly-Williams method within the tetrahedron approximation. The dependence 
of the results on the choice of different sets of ordered structures is discussed. 
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Introduction 


Recently, a new semiempirical method for calculating defect energetics in metallic alloys 
was introduced by Bozzolo, Ferrante and Smith (BFS) 1 . This technique, which builds on 
the ideas of Equivalent Crystal Theory (ECT) 2 was successfully applied to the study of heat 
of formation and lattice parameters of fee alloys as a function of alloy composition. BFS is a 
quantitatively accurate and computationally simple technique for determining the energetics 
of ordered multicomponent structures. Although there has been extensive calculations for 
fee alloys, similar results for bcc alloys have been limited. This is in part due to possible 
limitations in application of such approaches to bcc metals. Since alloys of bcc metals are 
important in structural materials, the present work represents an important contribution 
to the calculation of defect energies in alloys. In this effort, we apply BFS to the study of 
bcc-based binary alloys using the method of Connolly and Williams (CWM) 3 for the study 
of the energetics of disordered structures within the tetrahedron approximation. 

In Section 2 we present the BFS method and discuss the application of CWM to several 
choices of ordered structures. An application to selected bcc-based binary alloys is discussed 
in Section 3. Conclusions are drawn in Section 4. 
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Formalism 


In BFS, the energetics of binary alloys is described in terms of pure metal properties and 
only two experimentally (or theoretically) determined alloys properties. We build on the 
formulation of ECT by dividing the total energy of the alloy into a chemical energy and a 
strain or structural energy. The strain energy associated with a given atom is computed 
as if all of its neighbors were of the same atomic species. It arises from neighbor locations 
being different from in the elemental single-crystal environment. The remainder of the total 
energy is defined to be the chemical energy, which is due to some of an atom’s neighbors 
being of a different atomic species. We now proceed to outline the procedure for calculation 
of heats of formation versus concentration for alloys with multiple atomic species. With 
this procedure, the binding energy curve as a function of volume is obtained from which 
the bulk properties of specified alloys can be extracted. The application of this technique 
to different crystallographic structures is straightforward. In this work, we concentrate on 
bcc-based binary alloys. 

Consider a cell containing Nx atoms of atomic species X, ( X = A y B y ...), so that the 
total number of atoms in this cell is given by N = Yx Nx* The heat of formation of this 
cell is 

A E ce u = E ce ll - ^ ~2,NxEx (1) 

X 

where E ce u is the total energy of the cell and Ex is the cohesive energy of an atom of species 
X in a pure crystal of its own species. If E'(i,X) denotes the energy of the i-th atom in 
the cell (i = 1, Nx) of species X then 

N X 

= £!>,*. ( 2 ) 

X i= 1 
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where the energy difference €{ x — E\i, X) - Ex, has a strain and a chemical energy contri- 
bution, linked by a coupling factor g± t x that ensures that the chemical energy contribution 
vanishes for large interatomic distances: 

e i,x = ef,x + 9i,xe? tX ( 3 ) 


In order to compute the strain energy, we just ’flip’ every atom surrounding atom 
{i,X) into an atom of the same species X, and perform a regular ECT calculation 2 . The 
equivalent lattice parameter af x is determined by solving the appropriate ECT equation 
applied to atom (i, X) in the defect (but pure) crystal. Then, 


4x = E$[ 1 - (1 + o?jt)e-fr] , afx = (af x - a?)/l x . 


(U 


where a * denotes the lattice constant of a pure X crystal and Eq the corresponding cohesive 
energy. 

In Eq. (3), the strain and chemical energies are coupled nonlinearly. The coupling func- 
tion gi t x guarantees that the chemical contribution will vanish with increase in interatomic 
spacings, as it should. We define the coupling function in terms of the scaled equivalent 
lattice parameter of the strained crystal as follows 


9i,x = e 


5 * 

i,X 


( 5 ) 


For the chemical energy contribution we keep the actual chemical composition 

of the cell (i.e., A, B, C in proper proportions for the alloy) , but we force the atoms 
surrounding atom (i,X) to be located in the lattice sites of a pure crystal of species X. 
Thus, we are including the effect of changing a neighbor to an A, B or C atom). Two similar 
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ECT calculations are then performed: 


e ?,X = e< i,x{{&Y,x}) ~ e&c(0). (6) 

For the first term, e^({Ay t x)}, the chemical perturbation is included in the appropriate 
values of the set of parameters {Ayj} which include the effects of changing the atomic 
species of a neighbor, where Y denotes the atomic species of a given neighbor of atom 
(i,X). In order to determine the parameters A y,x and A x,Y for a given pair of atomic 
species (X,Y), two experimental values of any property of the X x Yi~ x alloy are needed. 
We choose to use the experimental heats of solution in the dilute limit which in most cases 
are readily available. 

The equivalent lattice parameter af x is obtained by solving the corresponding ECT 
equation 2 

12 R p 1 x e~ axRi +QR p 2 x e ~ {ax+ ^ )R2 = £ NxYr p 1 x e-^ x+AYX ') rt 

Y 

+ ^2M XY rYe~ {ax+ ^ +AYx)r2 ( 7 ) 

Y 

where R\ and R2 are the nearest-neighbor and next-nearest-neighbor distances in the equiv- 
alent crystal of lattice parameter af x . The first term in the chemical energy, ep^(A y,x) is 
then given by 

<x(Ay,*) = 7 i,xE$ [l - (1 + a?; x )e<'*\ , a& = (og* - of )fl x (8) 

with 7 i >x = 1 if af % > 0 and 7 j fX = — 1 otherwise. The second term in Eq(8) is obtained 
by a similar procedure, but setting all the perturbative parameters {Ay-.x} equal to zero. 
This is done in order to free the chemical energy from any structural defect information, 
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retaining only the contribution of the chemical composition of the surroundings of atom 

(«*.*)• 


In this work, we are concerned with calculating the heats of formation of the bcc-based 
ordered binary alloys A X B\_ X . If no relaxation of the individual atomic sites is allowed, 
then the strain energy is simply 

e x = [l ~ (1 + a x*)] e ~ a * » a x* = ( r “ a f )/(*• (9) 

where r is the actual interatomic distance of the alloy. Within this approximation, the 
second term in the chemical energy (Eq.(6)) vanishes, leaving us only with the computation 
of the first term, e^(A). For a given ordered structure m, the ECT equation for the 
equivalent lattice parameter a (Eq.(8)) is 

12 R\ x e~ axRl + 6 Rl x e' {ax+ ^ )R2 = Nxxr^ x e~ axri + N X yr^ x €~^ x+AYX '> n 

+ Mxxr?e-( ax+ ^ )r > 

+ Mx Y r p 2 x e- {ax+ ^ +£kYx)ri ( 10 ) 

with R\ = &R 2 ;R 2 = = ^ r 2> r 2 = a e • The parameters px , ax and \x are listed 

in Ref. 2 and the coefficients Nxx^XY^XX and Mxy depend on the different ordered 
structures considered. 

As in previous applications, we will use the heats of solution in the dilute limit as the 
experimented input for determining the parameters Aab and Aba • In order to compute the 
heats of formation of the disordered alloys, we use the Connolly- Williams method 3 . This 
method is based on a formal expression for the total energy first derived by Sanchez 4 , where 
the total energy of a particular configuration m of a binary alloy consisting of atoms A and 
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B on a lattice of fixed symmetry is given by 


AiUM = (11) 

7 

where v 7 (r) are many-body potentials, the are multisite correlation functions defined on 
a 7-type cluster, r is the lattice parameter and the sum includes all 7-type clusters on the 
lattice. The multisite correlation functions are defined as 

= 122 Vru-'-Pny ( 12 ) 

7 {«.} 

where <r n is a spinlike variable which takes the values +1 and —1 depending on whether 
the lattice point n is occupied by an A or B atom , and N y is the total number of 7-type 
clusters. 

The many-body potentials v y (r) are obtained by inversion of Eq. (11), which implies 
the existence of a maximum cluster 'y max beyond which the n 7 (r) are supposed to be neg- 
ligible. Thus, for a certain set of ordered structures a and by arbitrarily truncating the 
summation in Eq. (12), the many-body potentials are 

v ~/( r ) = ^ (c) & E m(r) , <f> < 7 < 7 max 

m 

V^,(r) = 0 , Tmax ^ 7 ^ 00 (13) 

where represents the empty cluster. Recently 5 , the CWM was extended to include more 
ordered structures and cluster sizes than the ones originally proposed 3 . Multisite correla- 
tions for the most common bcc and fee based superstructures were also given 5 . Table 1 lists 
the correlations included in the tetrahedron truncation of the CWM for some structures 
on the bcc lattice. Table 2 shows the coefficients Nxx > Mxx and Mxy needed to 

solve Eq. (10) for all the possible ordered structures included in Table 1. These ordered 
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structures are derived from the tetrahedron cluster shown in Fig. 1: a and 7 are on body 
centers and 0 and 6 are on cube edges. When a — 7 and 0 = 6 the structure is called B 2. 
The B32 structure is derived when a = 6 / 0 = 7 and the DO3 structure is obtained when 
(3 = S / a ^ 7. 

In this work, we considered different choices of ordered structures, as well as the type of 
clusters included in Eq. (12). Being that the experimental input is, obviously, the same for 
all cases studied, a comparison with available experimental data for the heats of formation 
of binary alloys 6 should give us an indication of the preferred ordered structures for a given 
binary alloy. 

The different choices are related to the two possible ordered structures at 50 % compo- 
sition ( B2 and B32) and the corresponding pair multisite correlation functions (£2 and £3). 
We will denote the cases studied as follows: (i) B32/^2* includes the B32 ordered structure 
and the £2 correlation function (nearest-neighbor pair), leaving out the B 2 structure as well 
as the £3 function (next-nearest-neighbor pair); (») B32/£3: includes the ordered structure 
B32 and the £3 function, leaving out B2 and £2; (in) B2/£2: includes B2 and £ 2 , leaving 
out B32 and £3 and (iu) B2+B32/£2 + £3: includes all the structures and functions listed in 
Table L In each case, the excess energy A E m (r) for the corresponding ordered structures 
is obtained with Eq. (2). Within the tetrahedron approximation, this calculation involves 
just a few atoms, as indicated in Table 2. 

Following CWM, the excess energy for the disordered alloys A X B\_ X is given by 

AE D (r,x) = ]jT(l - 2ar) n ''u Y (r) (14) 

7 

where n 7 is the number of sites contained in the 7 cluster. For each choice of ordered 
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structures and multisite correlation fuctions we have different many-body potentials 

(r) = J2(^y lAE m(r)- (15) 

m 

Replacing Eq. (15) in Eq. (14), we can write A Eo(r,x) for each one of the cases studied 
as 

AE D (r, x) = ^2 c m (x)AE m (r) (16) 

m 

where the sum runs over the appropriate ordered structures included in each case considered 
and the polynomials c m (x ) are also dependent on the clusters and structures included in 
each case. Table 3 lists the polynomials c m (x) for the reduced basis sets (t), («) and (fit) 
and Table 4 displays the corresponding polynomials for the general case (all structures and 
multisite correlation functions included in Table 1). Finally, the heat of formation for a 
given concentration x is obtained by finding the minimum value of AEo(r, x). 
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Results and Discussion 


In this section we present results for selected bcc-based alloys which display quite different 
behavior. For the four systems studied, we used the experimental values of the heats of 
solution in the dilute limit 6 , listed in Table 5. Table 6 displays the values of p, a, / and Ec 
for the pure elements 2 , needed to solve Eqs. (4)-(10). Table 7 shows the values of Aab and 
Aba one obtains with our procedure for the different choices of basis sets described in the 
previous section. 

The parameters A ab and Aba can be taken as ’perturbations’ to the pure-element 
a’s listed in Table 6, trying to simulate the interaction between two atoms of different 
species. In all cases, these quantities are small compared to the pure-element a’s, and, 
surprisingly, rather insensitive to the different choices of basis sets. However, these small 
differences translate into a noticeable change in the heat of formation versus concentration 
curves obtained for each choice of basis set, as can be seen in Fig. 2. 

Fig. 2. a shows the results obtained for Cr-Fe alloys where the regular, symmetric 
behavior of the heat of formation curve is accurately reproduced by using the basis set (»*»), 
where a B2 structure is included. Although there is no known ordered phase of Cr-Fe with 
the B2 structure, the fact that the choice (in) is favored over the others can be taken as 
an indication of the possible existence of such phases. This assumption is further validated 
by similar results for Fe-V alloys (Fig. 2.b) where, again, basis set (Hi) best approximates 
the available experimented values of the heat of formation 6 . In this case, there is some 
experimental evidence that such an ordered phase exists 6,7 . 

The predictions for Cr-Mo alloys show a drastic change in behavior as compared to 
Cr-Fe, although both systems display similar features in the experimental values of the heat 
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of formation. Fig. 2.c shows the theoretical results. In this case, the basis sets ( i ) and (ii) 
yield comparable results, predicting a preference for a B32 structure, although there are no 
known ordered phases of these alloys. 

Of the four examples shown in this work, Cr-V (Fig. 2.d) displays the most surprising 
features, therefore providing a severe test to the sensitivity and accuracy of our method. In 
the large body of experimental data for binary alloys, Cr-V is one of the very few to display 
the irregularities seen in the heat of formation vs. concentration curve, characterized by a 
sudden change in curvature for a small range of concentrations. As it was the case for Cr- 
Mo, this feature of Cr-V is approximately reproduced only by the results generated with the 
choice of the basis set (it). Once again, no ordered phase is known, but the clear distinction 
between the different basis set choices shown in Fig. 2.d can be taken as an indication of 
the crystallographic structure of these alloys. 

The fact that we used the experimental heats of solution in the dilute limit (i.e., the 
derivative of the heat of formation at x = 0 and x = 1) might lead one to believe that that 
choice somehow predetermines the behavior of the heat of formation curves. The examples 
shown in this work obviously contradict this fact: all four curves (for each choice of basis 
set) were obtained with the same values of the heats of solution. However, their behavior 
for the whole range of concentration is quite different in each case. 

The explanation for the particular features of the heat of formation vs. concentration 
curves is not then in the heats of solution, which is our only experimental input, but in 
the delicate balance between the strain and chemical energies, as defined in our formalism. 
Except for the case of Cr-Fe, where the small lattice mismatch results in an almost negligible 
strain energy, in all the other cases the heat of formation predicted with our method is 
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obtained from large competing strain and chemical energy contributions. 

Fig. 3 displays these contributions for the four systems considered in Fig. 2, showing 
the results obtained with the best basis set choice for each system, as discussed before. The 
apparent similarity seen in Cr-Fe and Cr-Mo for their heats of formation, arises from quite 
different strain and chemical energy contributions: while the positive chemical energy is 
mainly responsible for the heat of formation for Cr-Fe alloys (Fig. 3. a), a large negative 
chemical energy in Cr-Mo (Fig. 3.c) is necessary to balance a large strain energy contribu- 
tion, absent in Cr-Fe. Also, the symmetry seen in the heat of formation curve for Cr-Mo 
is a result of completely different regimes in the strain and chemical energies: linear for 
Mo-rich systems and with a pronounced curvature for Cr-rich alloys. Fig. 3.b shows, for 
Fe-V, how the chemical energy is solely responsible for the axis-crossing seen in the heat of 
formation curve. 

As noted before, Cr-V provides the appropriate grounds for testing the sensitivity 
of this method. Fig. 3.d displays the strain and chemical energy contributions for the 
Cr-V systems. One can see how a barely noticeable flattening in the chemical energy 
contribution is the source of the unusual feature seen in the heat of formation curve. These 
results correspond to the basis set (ri), which best approximates the experimental results. 
Finally, Fig. 4 expands on the results shown in Fig. 3.d in that the strain and chemical 
energy contributions are displayed for all four basis sets. The strain energy contribution 
(independent of A CrV and A vcr) shows small differences due to the choice of different 
ordered structures. The chemical energy term dictates the behavior of the heat of formation 
as a function of concentration: the asymmetry seen in Fig. 3.d arises from the chemical 
behavior of the B32 structure. 
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Conclusions 


In this work, we applied the semiempirical method of Bozzolo, Ferrante and Smith to the 
study of bcc-based alloys. The method was used to compute the total energy of ordered 
structures. The energetics of disordered alloys was studied with the CWM and several 
choices of basis sets were considered. Good agreement with experimental results for certain 
choices of basis sets gives an indication of the possible symmetries underlying the ordered 
compounds. The partition of the heat of formation into strain and chemical contributions 
provides some insight in the physical behavior of the systems studied. 
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TABLE 1. - CORRELATIONS INCLUDED IN THE TETRA- 
HEDRON TRUNCATION OF THE CWM FOR SOME 
STRUCTURES ON THE BCC LATTICE 

[The £j with increasing index i, correspond 
to the empty cluster, the point, the nearest- 
neighbor (NN) pair, the next-nearest-neighbor 
(NNN) pair, the triangle formed by two NN 
pairs and one NNN pair, and the tetrahedron 
formed by four NN pairs and two NNN pairs.] 


Structure 

Composition 

to 

Cl 

£2 

& 

u 

£5 

bcc 

A 


1 

1 

1 

1 

1 

do 3 

A 3 B 



0 

0 

1 

-1 

B2 

AB 



-1 

1 

0 

1 

B32 

AB 



0 

-1 

0 

1 

do 3 

ab 3 



0 

0 

1 

2 

-1 

bcc 

B 



1 

1 

-1 

1 


TABLE 2. - COEFFICIENTS N xx , N XY , M xx , and M XY 


Struc./Comp. 

Atom 

N>i/i 

N AB 

Nbb 

Nb.4 


M>jb 

Mbs 

Mb4 

bcc /B 



0 

8 



0 

6 

0 



mm 

8 

0 

0 

0 

6 

0 

0 

t>o 3 /ab 3 

B(2) 

mm 

0 

4 

4 

0 

0 

6 

0 


B(l) 

wM 

0 

8 

0 

0 

0 

0 

6 

B2 /AB 

A(2) 

0 

8 

0 

0 


0 

0 

0 


B(2) 

0 

0 

0 

8 


0 

6 

0 

B32 /AB 

EBB 

4 

4 

0 

0 

0 

6 

0 

0 


BBl 

0 

0 

4 

4 

0 

0 

0 

6 


wsm 

4 

4 

0 

0 

6 

0 

0 

0 

DO 3 /A 3 B 

m I 

8 

0 

0 

0 

0 

6 

0 

0 


189 

0 

0 

0 

8 

0 

0 

0 

6 

bcc /A 

A(4) 

8 

0 

0 

0 

6 

0 

0 

0 


TABLE 3. - POLYNOMIALS c m (x) FOE CASES (i), (ii), AND (iii) 


Struct. 

Comp. 

0 

(0 

(») 

(iii) 

bcc 

B 

0 

1 — Ax + 5x 2 — 2x 3 

1 — 4a; + 7a: 2 — 6x 3 + 2a: 4 

1 — 4x + 7x 2 - 6x 3 + 2x 4 

do 3 

ab 3 

1 

4x — \2x 2 + 12x 3 — 4x 4 

4x — 12x 2 + 12a: 3 - 4x 4 

4x - 12x 2 + 12x 3 - 4x 4 

B2 

AB 

2 

- 

- 

4x 2 — 8x 3 + 4x 4 

B32 

AB 

2 

8x 2 - 16a: 3 + 8a: 4 

4x 2 — 8x 3 + 4x 4 

- 

do 3 

A 3 B 

3 

4a: 3 - 4a: 4 

4x 3 - 4x 4 

4x 3 - 4x 4 

bcc 

A 

4 

-x 2 + 2a: 3 

x 2 - 2x 3 + 2x 4 

x 2 — 2x 3 + 2x 4 
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TABLE 4. - POLYNOMIALS c in (x) FOR CASE (iv) 


Structure 

Composition 

m 

(iv) 

bcc 

B 

0 

1 — 4x + 6a: 2 - 4x 3 + x 4 

do 3 

ab 3 

1 

4x — 12x 2 + 12a; 3 — 4a: 4 

B2 

AB 

2 

4a: 2 — 8a: 3 + 4a: 4 

B32 

AB 

2 

2a: 2 - 4a: 3 + 2a: 4 

do 3 

a 3 b 

3 

4x 3 — 4x 4 

bcc 

A 

4 

x 4 


TABLE 5. - EXPERIMENTAL 
HEATS OF SOLUTION 
E An AND E ra 


TABLE 6. - INPUT PARAMETERS FOR BCC ELEMENTS 


Element 


Cohesive 




Lattice 
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Figure 1.— Tetrahedron cluster in a bcc lattice (see text). 
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(a) Cr-Fe. 


(c) Cr-Mo. 


Figure 2. — Heat of formation versus concentration for different bcc-based alloys. In all cases, the solid squares indicate 
experimental values, and the different curves the results obtained using the basis sets and correlation functions 
described in the text (i) B32 1$ 2 (long-dashed line), (ii) B32/£ 3 (short-dashed line), (iii) B2/£ 2 (solid line) and (iv) B32+B2/£ 3 
(dotted line). 



















Chemical Strain 



Figure 4. — Strain and chemical contributions to the heat of formation 
of Cr-V for different choices of ordered structures (see text): 

(i) (longdashed curves), (ii) solid curves, (iii) (short-dashed curves 
andfiv) (dottedcurves). The sum of these pairs of curves are shown 
in Fig. 2.d. 
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